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A SIMPLIFIED METHOD FOR THE DETERMINATION AND ANALYSIS OF 
THE NEUTRAL-LATERAL-OSCILLATORY-STABILITY BOUNDARY 

By Leonard Sternfield and Ordway B. Gates, Jr. 


SUMMARY 

A necessary condition for neutral oscillatory stability is that 
Routh’s discriminant R, formed from the coefficients of the 
stability equation, is equal to zero. The expression for R is 
D(BC — AD) — B 2 E where A, B, C, D, and E are the coefficients 
of the lateral-stability equation. In a large number of the 
cases considered in this study, it has been found that the term 
B 2 E may be neglected. Routh’s discriminant is then factorable 
into two simplified expressions, that is, BC—AD=Ri and D; 
and either R\~0 or D=0, or both, may constitute a condition 
of neutral stability. Test functions have been derived which , if 
satisfied, indicate, that the simplified expressions may be used to 
approximate R=0. If R, =0 and D—0 satisfy the necessary 
and sufficient conditions for a neutral-oscillatory-stability 
boundary, D—0 represents the boundary for the oscillation 
which has a period comparatively longer than the period, of 
oscillation for which Rx—0 is the boundary. 

In general, the results of the computations obtained from 
Ri=0 and D=0 show very good agreement with the results 
calculated by the exact expression for 7?=0. The nature of the 
modes of motion as a function of the directional-stability deriv- 
ative and the effective-dihedral derivative is discussed in detail. 

INTRODUCTION 

The results of recent investigations (references 1 and 2 
and unpublished results of lateral-stability analyses for 
several experimental high-speed airplanes) have indicated 
that small variations in some of the airplane mass and 
aerodynamic parameters may cause a pronounced change in 
the oscillatory stability of the airplane. It has been difficult 
to explain the reasons for such pronounced changes because of 
the complexity of the expression for neutral oscillatory stabil- 
ity. This expression, based on the lateral-stability equations 
with three degrees of freedom, involves a large number of 
combinations of the mass and aerodynamic parameters. In 
order to predict the stability of the lateral oscillation, there- 
fore, it appears necessary to make a separate stability 
analysis for each airplane. 

The simplified expressions derived for the neutral- 
oscillatory-stability boundary in the present theoretical inves- 
tigation simplify the calculations required to obtain the bound- 
ary in the analysis essential for each airplane. Because of the 
comparative simplicity of these expressions, an insight into 

866769—50 


the important combinations of mass and aerodynamic 
parameters that affect the lateral oscillatory stability is also 
provided. Through further investigation and analysis of 
the effects of these major parameters, the necessity of making 
separate calculations for each airplane might possibly be 
eliminated. Test functions are given which, if satisfied, 
indicate that the simplified expressions may be used. 

The nature of the modes of motion as a function of (?„ 
and Ci ff , the directional-stability derivative and effective- 
dihedral derivative, respectively, is shown to depend upon 
the location of the stability boundaries plotted as a function 
of C nfj and C tp . 

The results of the calculations based on the simplified 
expressions are presented for comparison with the results 
obtained by the complete expression for the neutral- 
oscillatory-stability boundary. 

SYMBOLS AND COEFFICIENTS 

4> angle of bank, radians 

\p angle of azimuth, radians 

8 angle of sideslip, radians (v/V) 

v sideslip velocity along the F-axis, feet per second 

V airspeed, feet per second 

p mass density of air, slugs per cubic foot 

q dynamic pressure, pounds per square foot ^ pV 

b wing span, feet 

S wing area, square feet 

W weight of airplane, pounds 

m mass of airplane, slugs (W/g) 

g acceleration due to gravity, feet per second per second 

Pi, relative-density factor ( m/pSb ) 

r\ inclination of principal longitudinal axis of airplane 

with respect to flight path, positive when principal 
axis is above flight path at the nose, degrees 
(see fig. 1) 

9 angle between reference axis and horizontal axis, 

positive when reference axis is above horizontal 
axis, degrees (see fig. 1) 

€ angle between reference axis and principal axis, 

positive when reference axis is above principal 
axis, degrees (see fig. 1) 
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K 

K Xo 

Kz 0 


K x 


K z 

K xz 

C L 


C \ 


Cn 

Cy 

Ci* 


Cy ■„ 




c, 


V 


c h 


Cy, 


Cy 


t 


D„ 

R 


angle of flight path to horizontal axis, positive in a 
climb, degrees (see fig. 1) 

radius of gyration in roll about principal longitudinal 
axis, feet 

radius of gyration in yaw about principal vertical 
axis, feet 

nondimension al radius of gyration in roll about 
principal longitudinal axis (k x Jb) 
nondimensional radius of gyration m yaw about 
principal vertical axis ( k z Jb ) 
nondimensional radius of gyration in roll about 
longitudinal stability axis (yK x * cos 2 i>+ J K z 2 sin 2 ij) 
nondimensional radius of gyration in yaw about 
vertical stability axis (y'Kz 0 2 cos 2 t]+Kx 0 2 sin 2 rj) 
nondimensional product-of-inertia parameter 
({K z 2 —It x o 2 ) sin v cos n) 


( C OS y\ 

\qS-y 


trim lift coefficient 
rolling-moment coefficient 
yawing-moment coefficient 
lateral-force coefficient ( 


/■Rolling moment^ 

V qSb ) 


qSb 

/ Yawing moment\ 

v' " qSb ' ) 

/Lateral force\ 
v qS / 
effective-dihedral derivative, rate of change of rolling- 
moment coefficient with angle of sideslip, per 
radian (dC;/d/3) 

directional-stability derivative, rate of change of 
yawing-moment coefficient with angle of sideslip, 
per radian {dC n /d0) 

lateral-force dei'ivative, rate of change of lateral- 
force coefficient with angle of sideslip, per radian 
(pCy/dfi 

damping-in-yaw derivative, rate of change of yawing- 
moment coefficient with yawing-angular-velocitv 

factor, per radian ~ j 

rate of change of yawing-moment coefficient with 
rolling-angular-velocity factor, per radian 

(acyaf/) 

damping-in- roll derivative, rate of change of rolling- 
moment coefficient with rolling-angular-velocity 


factor, per radian 


(d6',/d |/ 


rate of change of rolling-moment coefficient with 
yawing-angular-velocitv factor, per radian 

( bC ^Tv) 

rate of change of lateral-force coefficient with rolling- 
angular-velocity factor, p or radian (bC Y [b 
rate of change of lateral-force coefficient with yawing- 
angular-velocity factor, per radian (i>C Y ! d 
time, seconds 

nondimensional time parameter based on span (Vt/b) 
differential operator 
Routh’s discriminant 


\ 



Fkh'rb I. — System of axes and angular relationships in flight. Arrows indicate positive 
direction of angles. V=0—y—*. 

X complex root of stability equation 

AX 4 + B\ 3 + C\ 2 -\-D\-\-E=0 (X— £drico) 

X' complex root of stability equation 

2x , 4 +Rx , 3 +cx , 2 +I>x , +£ i =0 (x'=r±f«0 

P period of oscillation, seconds 

T U2 time for amplitude of oscillation to change by factor 
of 2 (positive value indicates a decrease to half- 
amplitude, negative value indicates an increase to 
double amplitude) 

A, R, C, D , E coefficients of lateral-stability equation 
EQUATIONS OF MOTION 

The nondimensional linearized equations of motion, re- 
ferred to the stability axes, used to calculate the spiral- 
stability and oscillatory-stability boundaries for any flight 
condition arc: 

Rolling 

2» b {K x *DU+K xz DW) C lp D^+ \ C t D ^ 

Yawing 

‘2^b(K z 2 D 2 \l' J rK X zDi, 2 4>) = Cngfi-i-^ C „ pD C V t D b 4> 
Sideslipping 

2fi„(D b l3 J rI) b il/)= C Y jL) b <b+C L <t) J r-^C Yr D b 4/ J r{CL tan 7)/ 

When </> 0 e XSh is substituted for <j>, ft 0 e Ub for 0 o e XSk for 0 in 
the equations written in determinant form, X must be a root 
of the stability equation 

.4X 4 +RX 3 -(-(7X 2 +Z>X+E'=0 (1) 


where 

A=^{K x 2 K-/-K X7 2 ) 

B= — 2^(2K x -K/C V0 +K x 2 C nr + K z HJ lp - 2 K xt *C Yft - 
K xz C lr -K xx C np ) 


J 
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c= p b ( K x *C nr C Yf) + Ap b K x 2 C ng + KzC \ v C Y p + \ C ri C lp - 

KxzClCyp 4:p b K XZ C C n K x zC Yff — ^ n lP l T 

K X z G H p Gy p Kz 2 Cy p C Ip — K X 2 Cyfinp^ K X Z Gyfii^J 

D— — 4 Cn r C lp Cyp — p t Gl p Gnp + -^ Cn p Gi r Cyp+p l ,G„ ii C,p + 
2ptCi,KxzCnp 2p 0 C L K z 2 Ci fj — 2p b K x 2 GnpC L tan y-\- 
2p b KxzC i ff C L tan 7 +^ Ci p C n pCy r — -^ G , lp G i p G Y r — 
\C lr C nB Gy v +\Cn T G l pGy v 

E=\ C L {C nT C lB -C lr C n p)+\ C L tan pC.C^-C.C^) 


Multiplying equation (1) by p b and substituting X= - 
results in the stability equation 


where. 


A\ n + B\' 3 + C\ ri + D\' + E= 0 


7 A 


A ii 

M& 


( 2 ) 


5=4 

Pb 


C— - 

Pb 


D=D 


E=p„E ■ 

The damping and period of the lateral oscillation in 
seconds are given respectively by the equations 

rp __ 0.69 p b b 

-ti/2 £/ y 

p 6.28 p„b 

r u' V 


where tj' and to' are the real and imaginary parts of the 
complex root of stability equation ( 2 ). 

ANALYSIS 


If X=fw is substituted in the stability equation, the following 
two equations are obtained: 


Au 4 -Gco 2 +E=0 (4) 

-Ba> s +Do> = 0 (5) 


Solving equation (5) for o > 2 and then substituting the result 
("'“§) in equation (4) results in Routh’s discriminant 
B CD —AD 2 — B 2 E— 0 


It is seen from equation (5) that 


VI 


is the angular fre- 


quency of the neutrally stable oscillation only when B and 
D are of the same sign since cc must have a real value if the 
root X=±ho is to represent an oscillation. If B and D are 
of opposite sign and R = 0 is satisfied, the two roots of the 
stability equation given by X=±iw are real, one positive 
and one negative. It is important to note that the A, G, 
and E coefficients may be of opposite sign to the B and D 
coefficients, and neutral oscillatory stability will still occur 
as long as Roulb’s discriminant is equal to zero and the D 
and B coefficients are of the same sign. 


In genera], the R= 0 boundary in the G n/} ,Gi g plane lias 
two branches. The two branches result from the fact that 
7?=0 is a quadratic equation in C, ff and thus lias two C t0 
roots for every value of C np . Usually, the two branches can 
be approximated by simplified expressions for R= 0. In 
certain cases, however, which arc discussed in the section 
entitled “Test Functions,” either one or none of the branches 
may be approximated. 

Now, the condition f?=0 is a necessary but insufficient 
condition for neutral oscillatory stability. The simplified 
expressions, therefore, which approximate R — 0 do not 
necessarily represent boundaries of neutral oscillatory 
stability. Other conditions, elaborated on in the section 
“Validity of D= 0 and Ry = 0 as Ncutral-Oscillatory-Stability 
Boundaries,” must bo satisfied before either expression 
represents a valid boundary. 

There are, therefore, two kinds of tests that must be made: 
First, tests to determine whether R= 0 may be approximated 
by simplified expressions; and, second, tests to determine 
which of these expressions represents a boundary of neutral 
oscillatory stability. The significance of the lateral-stability 
boundaries is indicated by a discussion of the modes of motion 
in the G„ B ,Gi p plane. 


The necessary and sufficient conditions for neutral oscilla- 
tory stability, as shown in reference 3 , are that the coeffi- 
cients of the stability equation satisfy Routh’s discriminant 
set equal to zero 

R = BCD AD 2 — B 2 E=0 ( 3 ) 

and that B and D have the same sign. The expression for 
R = 0 can be derived by assuming that the stability equation 
lias two roots X ~ ~ re, where oj is the angular frequency of 
the neutrally stable oscillation. This assumption is based 
on the fact that for the condition of neutral oscillatory 
stability the real part of the complex root must be zero. 


DERIVATION OF SIMPLIFIED EXPRESSIONS 

Two of the most important stability derivatives affecting 
lateral oscillatory stability are the directional-stability 
derivative C np and the effective-dihedral derivative CV. 
The boundary for neutral oscillatory stability is usually 
plotted as a function of these two derivatives with C„ fi as 
the ordinate and C lfl as the abscissa. The method used to 
obtain the neutral oscillatory stability boundary is first to 
substitute the values of the mass and aerodynamic param- 
eters of a specific airplane in the coefficients of the stability 
equation while maintaining C„ p and C tp as variables and 
then to calculate the Ci p roots of equation (3) for several 
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values of C ng . For a given value of C ng , the expression for 
R=Q is a quadratic equation in C ig that is of the form 

Q’iCip-\-biCip-\-Ci—§ 

For a very large number of cases, it has been found that the 
last term of Routh’s discriminant B 2 E, which contributes 
only to bi and c,, produces a negligible change in the expres- 
sion 

&r Cip + b i Cip -f- Ci = 0 

If, therefore, the term B 2 E is neglected, equation (3) for 
Routh’s discriminant reduces to 

R=D(BC—AD)=0 (6a) 


TEST FUNCTIONS 

The approximate discriminants Ri—0 and D= 0 are based 
on the assumption that B 2 E can be neglected when Routh’s 
discriminant is set equal to zero. Thus, the simplified 
expressions for the neutral-oscillatory-stability boundary, 
Ri and D, should not be used if including the terms B 2 E 
causes an appreciable change in the roots of R= 0. In 
appendix A test functions are derived which indicate the 
incremental change in the roots of R l = 0 and D— 0 due to 
the terms B 2 E. If certain conditions placed upon these test 
functions are satisfied, then 7?i and D can be used to calculate 
the R= 0 boundary. 

If, at a given value of C ng , the root of 7?i = 0 is denoted 
by O lg —r, the approximate deviation of this root from a 
root of R= 0 is given by 


Thus BC—AD—0 and D = 0 are the approximate factors of 
the quadratic equation 

a\C i g -^biC i g -\- Ci=Q 


Ar= _ ( 7 ) 

Ar r 1 d 1 (d-r)+e 1 Uj 

If Ar is small, then l?i = 0 is a close approximation to one 
branch of R— 0. A suitable criterion for this approximation is 


The expression BC—AD is henceforth called R x . 

In order to simplify the expressions for Ri and D, the 
expected range of values of mass and aerodynamic param- 
eters for high-speed aircraft were substituted in the coeffi- 
cients of the stability equation to determine which terms 
could be omitted without appreciably affecting the values of 
Ri and D. The following simplified coefficients were obtained: 

A=^(K X 2 K Z 2 -K XZ 2 ) 

B=- 2^{K x 2 C nr +2K x 2 K z 2 C Vg + C,K./) 

0=1*1, ^ 4:fl t K x 2 Cnp — il*bKxzOi g -{--^ C n /\ \- 

C,K,‘C r ,- 1 <V\) 

D = -2C l K/) (Ci-2CJt xz ) C ng 

E=\ C L (C„O lg -C\C ng ) 


The expressions for R x and D thus become 
R\ = ( A\K X z — A 2 A 3 ) C; 0 +[ — Kxz&AiClA' C i p K xz ) K X ~(A\ 
\jQ ( A i- G h K **)-C lr C„ p y= 0 (6b) 

D=A 2 Ci g — ( Ci v — 2 C l Kxz) C ng = 0 (6c) 

Ai = K x 2 C nr +2K x 2 K/Cy g +C h K z 2 'j 
A 2 =Cn p —2C L K z 2 (6d) 

A 3 =K X 2 K Z 2 -K XZ 2 


and 

where 


The simplified expressions R x = 0 and D= 0, as presented, 
are applicable only to conditions of level flight or to condi- 
tions of small angles of glide or climb. Simplified expressions 
for conditions of large angles of glide or climb can be derived 
by a procedure similar to the one presented. 


or 


AH ^ 


r 

5 


|Ar| JsO.Ol 


vhichever is the larger. 

Similarly if a root of D = 0 is denoted by Ci g = d, the 
ipproximate deviation of this root from a root of R = 0 is 
pven by 




e, (e—d) 


If Ad is small, then D = 0 is a close approximation to one 
branch of R= 0. A suitable criterion for this approximation is 


|A d\Z 


d 

5 


or 

|dA|^0.01 


whichever is the larger. 

The expressions for r u d u e u r, d, and e for use in equations 
(7) and (8) are 

v 1 = 8i*b(AiK xz A 2 A 3 ) 

di=t*bA 2 

e i — 2fi b A l 2 C L C„ r 

0 

e ~ ( Gi p)b 2 E=0 

where A u A 2 , and A 2 are defined in equations (6d). 

The value of C ng to be used inthese test functions should be 
selected from the probable range of C„ g of the airplane for 
which the lateral-stability analysis is to be made. Thus, the 
approximation of Ri = 0 and Z) = 0 to .R = 0 is determined in that 
region of the C ng ,Ci g plane pertinent to a particular analysis. 
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VALIDITY OF D = 0 AND Ri=0 AS NEUTRAL-OSCILLATORY-STABILITY 
BOUNDARIES 

As mentioned previously, for R= 0 to be a boundary for 
neutral oscillatory stability, the coefficients B and D must 
be of the same sign. The three predominant terms of B con- 
tain the factors — C„ r , — C Yg , and —G lp , respectively. For 
positive damping in roll, C tp is negative; and for positive 
weathercock stability ( C np positive), C Yfj and C„ r are negative. 
Thus, B is positive in the usual case where there is weather- 
cock stability and damping in roll. Therefore, D must gen- 
erally be positive in order that R= 0 be a neutral-stability 
boundary. If the exact boundary R= 0 has been calculated, 
it is merely necessary to plot D= 0 and R= 0 and note whether 
R=0 is located on the side of D= 0 where D is positive. A 
primary purpose of the present report, however, is to obviate 
calculation of the exact boundary by the use of simplified 
boundaries together with test functions. A method to 
determine the sign of D from the results of the simplified 
expressions is therefore presented in the following paragraph. 

For a given value of C n (selected from the probable range 
of C n$ of the airplane for which the lateral-stability analysis 
is to be made) let d be a value of C, for which D= 0 and d' 
be a slightly different value for which 


R=R X D — BdE— 0 


The substitution of C tfi =d' gives 


Did') 


B 2 E(d') 
Ri id') 


The sign of D at the R= 0 boundary {C lp =d') is therefore 
determined by the signs of E and R x at d'. But since d 
differs little from d', the signs of E(d) and R x {d) will be the 
same as the signs of E(d') and R x (d'), respectively. Hence 
the sign of D at R= 0 is the same as the sign of E/R x at D= 0 
(fig. 2) ; that is, 


D{d') 


. E{d) 
~Bi id) 


If the signs of E and R x arc the same, I) is positive and repre- 
sents a neutral-oscillatory-stability boundary; if E and R x 
are of opposite sign, D is negative and then represents a 
boundary for which the roots of the stability equation arc 
equal and opposite in sign. 


The preceding analysis is applicable for the large majority 
of cases where E (d') and Ri(d') are of the same sign as E(d) 
and Ri(d), respectively. For these cases, the D curve is 
widely separated from the E and R x curves. If the D curve 
is close to either the E or R x curve, the signs of E and R x 
should be determined at C t0 =d'. However, a very good 
approximation to the value of d' can be obtained by adding 
to d the value of A d calculated in the previous section entitled 
“Test Functions.” Hence, the sign of D is determined from 
the signs of E and R x at C tg =d+Ad. 

If the value of Ci 0 at which R x =() is substituted in D and 
the resultant sign is positive, l?i=0 is a neutral-oscillatory- 
stability boundary. 

It is interesting to note that for some aircraft, the D= 0 
curve, which approximates one branch of the R= 0 curve, is 
a neutral-oscillatory-stability boundary over one section of 
the curve and a boundary for equal and opposite real roots 
over the remaining section. This division of the D=0 curve 
into two distinct parts is caused by a change in sign of the D 
coefficient at some point on the curve. If, as has been found 
in a large number of cases, R x is positive for all values of C„ 0 
and C tfj on the D= 0 curve, the sign of the D coefficient de- 
pends only on the sign of E at these points. As shown in 
figure 3, therefore, the point of intersection of the curves D = 0 
and E— 0 is the point of separation of the D= 0 curve into 
two characteristically different sections. For points on the 
hatched side of E=0, the E coefficient is negative and, there- 
fore, the dashed part of D — 0 is a boundary of equal and 
opposite real roots. Conversely, on the unhatched side of 
E= 0, the E coefficient is positive and the solid part of .0=0 
approximates a 'boundary of neutral oscillatory stability. 

For small positive or negative values of C n and negative 
damping in roll, it is possible for B to be negative. A similar 
analysis is applicable to this case where now I) must be nega- 
tive to satisfy the necessary condition that I? = 0 is a boundary 
of neutral oscillatory stability. 

In general, when the simplified expressions are used to 
obtain a neutral-oscillatory-stability boundary, the procedure 
to be used is as follows: 

(1) For a given value of C n0 , selected from the probable 
range of C n0 of the airplane for which the lateral-stability 
analysis is to be made, calculate r and d, the Ci 0 roots of 
R x = 0 and D= 0, respectively. 



Figure 2— Validity of D= 0 as a neutral-oscillatory-stability boundary. 



Figure 3.— Effect of the position of the E= 0 boundary on the validity of D~ 0 as an approx 
mate neutral-oscillatory-stability boundary. 
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(2) Determine the value of A r and A d by substituting the 
results into the test functions. 

(3) If the criterions for A r and Ad as set forth in appendix 
A are satisfied, consider 7?i = 0 and D — 0 close approximations 
to the R= 0 boundary. 

(4) In order to determine the validity of 7?i = 0 as a 
boundary of neutral oscillatory stability, substitute the 
given value of C n ' and C^—r into the D coefficient. (If the 

resulting sign is positive, 7?i = 0 approximates a branch of the 
neutral-oscillatory-stability boundary.) 

(5) In order to determine the validity of D= 0 as a 
boundary of neutral oscillatory stability, substitute the given 

E 

value of C n(l and C tg =d into ^ . (If the resulting sign 

is positive, D— 0 apimoximates a branch of the neutral- 
oscillatory-stability boundary; if the resulting sign is nega- 
tive, D= 0 approximates a boundary of equal and opposite 
real roots.) 

NATURE OF MODES OF MOTION IN THE C„ C, PLANE 

n F 'll 

In this section, the changes in the roots of the lateral- 
stability equation, which occur upon crossing the various 
stability boundaries, are discussed according to the principles 
of the theory of equations as given in references 3 and 4. 
The solution of the lateral-stability equation gives four roots 
which may be four real roots, two pairs of conjugate, complex 
roots, or two real roots and one conjugate complex pair. A 
pair of complex roots indicates an oscillatory mode and a 
real root indicates an aperiodic mode. If the airplane is 
disturbed from its trimmed condition by an arbitrary dis- 
turbance, the subsequent motion is compounded of these 
modes in different proportions. The method of calculating 
the different proportions of the modes is presented, for ex- 
ample, in references 5 and 6. Such calculations of the motion 
for numerous points throughout the plane would lie 

very laborious. It is more practical, therefore, to investigate 
merely the types of modes that may be expected throughout 
the C n .,C t plane as indicated by the stability boundaries. 

The calculation of the motion could then be limited to several 
points of interest. 

Consider a case where the neutral-oscillatory-stability 
boundary 7?i = 0 and the spiral-stability boundary E= 0 
are located in the first quadrant of figure 4 (a). The area 
between the two boundaries is a region of complete stability. 
The roots of the stability equation for combinations ol 
C ng and C tg in this region, such as point A in figure 4 (a), 
are two negative real roots and one conjugate complex pair 
with the real part negative. One of the real roots which is 
numerically small corresponds to the spirally stable motion 
of the airplane. The other real root corresponds to the 
heavy damping of the pure rolling motion. The complex 
roots with the real part negative show that the so-called 
Dutch roll oscillation is stable. Passing through the E= 0 
boundary from point A to point B causes the spiral mode 
to become unstable, and crossing through the i?! = 0 boundary 


from point A to point C causes the oscillatory mode to 
become unstable. The second branch of the R = 0 boundary 
plotted in the second quadrant as D — 0 is not a neutral- 
osciilatory-stability boundary but rather a boundary for 
equal and opposite roots as determined by the analysis 
presented in the section entitled “Validity of D = 0 and 
R 1 =0 as Neutral-Oscillatory-Stability Boundaries.” The 
roots of the stability equation for combinations of C n „ and 
C tg on this boundary are two equal and opposite real roots 
and a pair of complex roots with the real part negative. 
The positive real root is the spirally unstable mode, and the 
negative real root is the damping-in-roll mode. The 
oscillation continues to remain stable even though the D 
coefficient is negative. 

For the case where one oscillatory-stability boundary 
D— 0 appears in the first quadrant and another oscillatory 
stability boundary Ri — 0 is in the second quadrant (fig. 4 (b)), 
the period of the neutrally stable oscillation is much 
greater on D = 0 than on R { = 0. This fact can be shown to 
be true by investigating the angular frequency of the 
neutrally stable oscillation for points located on the i? 1 = 0 
and D = 0 boundaries. As shown previously, the angular 
frequency o> is equal to \L)jB; and, therefore, since the 
boundary I )= 0 approximates one branch of R= 0, the 
angular frequency for points on that branch is very small. 
For combinations of C n and C tfj on 7?i = 0, the angular 
frequency is much greater. In general, D= 0 is a neutral- 
oscillatory-stability boundary for a long-period oscillation. 
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The roots at point A of figure 4 (b) have the same character 
as the roots at point A of figure 4 (a), that is, two negative 
real roots and one pair of conjugate complex roots. At 
point B the roots of the lateral-stability equation are two 
pairs of conjugate complex roots. It is interesting to note 
that the boundary for two equal roots occurs between point A 
and point B and can be considered the boundary beyond 
which two pairs of complex roots exist. Reference 4 shows 
that for a quartic equation 

A\*+B\ 3 + C\ 2 +D\+ E= 0 


C ne and Cip where the complex roots divide into two real 
roots. 

There have been several cases for which a neutral-oscilla- 
tory-stability boundary did not exist in the C n ,C t plane. 
An analysis of these cases indicated that the boundary for 
equal roots was in the oscillatorily stable region and had 
divided the stable oscillation into two subsiding modes. 
The neutral-oscillatory-stability boundary, therefore, would 
not have any significance. 

RESULTS AND DISCUSSION 


the boundary for equal roots is obtained by setting the 
discriminant 

— 4P 3 — 27 Q 1 


equal to zero, where 

P = BD — AAE— ( '^ 

and 


Q=-B*E+ 


BCD 
3 + 


8 ACE 
3 


-AD 2 - 


2<7 3 

27 


Between this boundary and D— 0, the period of the stable 
oscillation which corresponds to the newly formed pair of 
complex roots is longer than the period of the. oscillation 
which corresponds to the other pair of complex roots. As 
Ci is increased to point C on the unstable side of D= 0, 
the newly formed long-period oscillation is the one that 
becomes unstable, whereas the short-period oscillation 
remains stable. At point D the roots consist of a spirally 
unstable mode, a stable mode due to the derivative (\ , and 

' *p } 

a stable oscillation which becomes unstable in passing through 
Ri — 0 to point E. 

Figure 4 (c) represents the case where both /?, = 0 and 
77=0 appear in the first quadrant but oidy R, is a neutral- 
oscillatory-stability boundary. The curve D= 0 is the 
boundary for two equal and opposite real roots. Point A 
once again has two real negative roots and a pair of complex 
roots with the real part negative. At point B, on the 
unstable side of 7?i = 0, the real part of the complex roots is 
positive and indicates an unstable oscillation, whereas the 
two real roots are still negative. The boundary for C— 0 is 
between 7?i = 0 and D= 0. Some investigators of lateral 
stability have thought that a radical change occurs in the 
roots upon crossing through this boundary. The calcula- 
tions indicate, however, that the roots do not vary appre- 
ciably upon passing through C= 0. At D= 0, however, 
there must exist two equal and opposite real roots; this 
condition is possible only if the complex roots divide into 
real roots since the other two real roots are negative in sign. 
The calculation of roots at point C indicate that the complex- 
roots had divided into two real positive roots, one of which 
was exactly equal in magnitude to one of the negative roots. 
Again, the boundary for two equal roots, located between 
<7=0 and point C, would determine the combinations of 


The simplified expressions were used to calculate i?i = 0 and 
D= 0, and the results are compared with the results of the 
calculation of R = 0 based on the complete expression. Not 
only do Ri — 0 and D = 0 show good agreement with 0 
(figs. 5 to 13) but the comparative simplicity of the Ri and 
D expressions allows identification of the major parameters 
that affect the stability boundaries. 

EFFECT OF C„ p -2Ci.Kz * ON THE BRANCH OF R = 0 
APPROXIMATED BYZ)=0 

Reference 2 shows that a stabilizing shift in the R = 0 
boundary is obtained when C, ly is increased in a positive 
direction up to a certain value, but further increases in the 
positive direction cause a destabilizing shift in R= 0. The 
effect of varying C„ on the R = 0 curve is presented in 
figure 5 for a model tested in the Langley free-flight tunnel. 
The figure, illustrates very good agreement between = 0 
and the simplified expressions J?i = 0 and Z7= 0. The 
expression for D = 0 is 

(Ci p 2 C l K X z ) C 

l0= ~ C , lp — 2 ClK-/ 

which indicates that for positive C n? when the numerator is 
negative in sign the D= 0 boundary is in the second quad- 
rant for negative values of C nji —2C^K Z 2 =A 2 and in the first 
quadrant for positive values of A 2 . For the cases of negative 
An presented in figure 5, the D = 0 boundary would appear 
in the second quadrant. It can be shown, however, by the 
method described in the section “Validity of Z7= 0 and Ri = 0 
as Neutral-Oscillatory-Stability Boundaries” that 77=0 in 
the second quadrant is not a neutral-oscillatory-stability 
boundary and hence is not plotted in figure 5. However, as 
C„ p is increased in a positive direction, where now A 2 is posi- 
tive, an increase in the positive value of A 2 causes the D = 0 
boundary to shift upward in the first quadrant in a destab- 
ilizing direction. 

From the results shown in figure 5, it is seen that for the 
cases of C np equal to 0.30 and 0.40 only the solid-line part of 
the R = 0 curve in the first quadrant (the branch which may 
be approximated by D= 0) is a neutral-oscillatory-stability 
boundary. The short-dash-line part of R=0 is a boundary 
of equal and opposite real roots. The reason for this division 
of the R = 0 curve into two parts is discussed in the section 
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Figure 5. — Effect of C» and -I 2 on the neutral-oscillatory-stability boundary. 


entitled “Validity of 77 = 0 and R l = () as Neutral-Oscillatory- 
Stability Boundaries” and is illustrated in figure 3. 

EFFECT OF C„ p -2CiKz 2 ON THE BRANCH OF R= 0 
APPROXIMATED BY Ri = 0 

The important effect of C„ p on 77! = 0 occurs only in the 
coefficient of the C lfj term, A X K XZ —A 2 A % , in which C„ p affects 
only the factor A 2 . The sign of A ?J is always positive and the 
sign of A\ is negative for positive G„ r By definition, K X z is 
positive if the principal axis is above the flight path at the 
nose of the airplane as is the case for the curves presented in 
figure 5. \In general, for positive C„ g , the expression of R x = 0 
which does not include any C, g terms is positive and, except 
for one term, is independent of C Hp . If, therefore, the 
coefficient of C !g is positive, 77i=0 is in the first quadrant; 
whereas if the coefficient of CV, is negative, R x = 0 is in the 
second quadrant. As C„ increases in a positive direction 
and A 2 becomes more positive, the coefficient of C’i g becomes 
more negative and Ri — 0 in the second quadrant shifts 
upward in a destabilizing direction. If A 2 is negative but 
the absolute value of A 2 increases, as in going from C„ p =0.15 
to C np =— 0.10 in figure 5, the coefficient of C lg becomes 
more positive and R\ in the first quadrant also shifts upward 
in a destabilizing direction. Thus the results indicate that 
increasing the absolute value of A 2 has a destabilizing effect 
on the neutral-oscillatory-stability boundary. 

According to a previous discussion herein, variations in 
C n that maintain A 2 constant cause no shift in the 77= 0 
boundary. When, therefore, 77=0 is approximated by 
77= 0, such changes in C n and K z should have a negligible 
effect on the 77= 0 boundary. In order to test this point, 


calculations were made for a free-flight airplane model for 
C n varying from 0.30 to 0.63 while simultaneously varying 
K z 2 in order to maintain the same positive value of A 2 . 
The results showed the expected insensitivity of the 77 = 0 
boundary to these changes. 

It should be remembered that 71=0 in the first quadrant 
is the neutral-oscillatory-stability boundary for the long- 
period oscillation; and if instability were to occur, the pilot 
might not find this type of instability difficult to control. 

EFFECT OF C„„ Cr g . AND Kx ON THE BRANCH OF R = 0 
APPROXIMATED BY D= 0 

The D expression indicates that the 77= 0 boundary is 
independent of the derivatives C, 1t and C Y e and the mas s 
parameter K x . Figures 6 and 7 show a comparison of the 
results obtained by the complete calculations with D= 0 
for the cases in which C Hr and C Yp , respectively, were arbi- 
trarily doubled in value. As noted in the figures, C„ r and 
C Yg have a negligible effect on the boundary. The effect of 
K x on the branch of 77=0 which may be approximated by 
D = 0 is shown in figure 8. Complete calculations were 
made to obtain the 77=0 curves for the previously discussed 
free-flight airplane model. The value of A 2 used in these 
calculations was 0.17. For purposes of comparison, K x was 
arbitrarily increased by a factor of 2.5. Again the results 
show practically no effect of K x on this branch of 77=0, as 
is indicated by the simplified expression 77=0. For the 
case discussed in figures 6 to 8, the branch of 77=0 approxi- 
mated by 77i = 0 is in the second quadrant and has little 
practical importance. Hence, the effect of these parameters 
on 77i was not determined for this particular case. 
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Figure 7 —Effect of Cy R on the neutral-oscillatory-stability boundary. 



Figure 8. — Effect of Kx on the neutral-oscillatory-stability boundary. 


EFFECT OF PRODUCT OF INERTIA ON THE BRANCHES OF #=0 
APPROXIMATED BY I2i=0 AND E= 0 

The product of inertia has been shown to have a very 
pronounced effect on the lateral stability of present-day 
airplanes designed for high-speed high-altitude flight (refer- 
ences 1 and 7). The importance of the product of inertia is 
illustrated in figure 9 (a), which presents the R — 0 boundaries 
of the hypothetical airplane discussed in reference 1 for two 
angles of inclination of the principal axis relative to the 
flight path, r?=0° and t]=2°. Calculations were also made 
for these cases using i?i=0; and the results presented in 
figure 9 (a) show the same marked stabilizing shift in the 
boundary, caused by the 2° inclination of the principal axis 
above the flight path, as obtained by the complete calcula- 
tions. The value of A 2 for the R i = 0 calculations was — 0.18. 
The value of C np was then increased so that A 2 was equal 
to 0.13 (fig. 9 (b) ) . In this case, D= 0 appears in the first 
quadrant and Ri = 0 is in the second quadrant. Although 
both D = 0 and R,=0 are valid boundaries, the results are 
discussed only for the effect of product of inertia on D = 0 
since only the C n ,C t combinations in the first quadrant 
are usually of practical significance. Calculations for ij = 0° 
and ij=2° were made using D = 0 and R= 0. Although the 
product-of-incrtia factor K X z does appear in the D expression 
(in the term —2C L K X z), an examination of D indicates that 
this term could have only a negligible effect on D = 0 when 



(a) Negative A 2 . 

(b) Positive A 2 . 

Figure 9.— E fleet of Kxz on the neutral-oscillatory-stability boundary. 
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Ci v is much greater than 2ClK xz , as is usually the case. 
Figure 9 (b) shows that the results predicted from D = 0 
agree very well with the results obtained from the complete 
calculations. 

EFFECT OF RADII OF GYRATION. ON THE BRANCH OF R=0 
APPROXIMATED BY Ri=0 

Figures 10 to 12 are presented for the purpose of showing 
the close agreement between results obtained by using Ri=0 
and results obtained from reference 1. The three figures 
illustrate the effect of the radii of gyration in roll and yaw 
k X(j and k z , respectively, on the neutral-oscillatory-stability 
boundary. Figure 12 emphasizes the fact that the simplified 
expression is sufficiently accurate to predict the effect of 
k Xa on the oscillatory -stability boundary throughout the 
entire range of variation of k Xn - 

EFFECT OF WING LOADING AND ALTITUDE ON THE BRANCHES OF R= 0 
APPROXIMATED BY R, = 0 AND D=0 

The effects of wing loading and altitude on the neutral- 
oscillatory -stability boundaries were determined simul- 
taneously by considering variations in the relative density 



-c,. 



Figure 11.— Effect of fc^ 0 on the* ncutral-oscillatory-stahility boundary. 



-Ci e 


Figure 12. — Effect of kx o on the neutral-oscillatory-stability boundary. 


factor because fi b varies directly with both wing loading 
and altitude. An examination of the expressions = 0 
and D = 0 indicated that increasing fi h causes a slight destab- 
ilizing shift in /?i = 0 but does not affect 79=0 since n„ does 
not appear in the expression for D= 0. The trend shown by 
these results agrees with the results found in reference 1 
concerning the effect of fi b on the neutral-oscillatory-stability 
boundary. 

COMPARISON BETWEEN NEUTRAL-OSCILLATORY-STABILITY BOUNDARIES 

OBTAINED BY EXACT AND SIMPLIFIED EXPRESSIONS FOR A HIGH-SPEED 

EXPERIMENTAL AIRPLANE 

Some of the neutral-oscillatory-stability boundaries ob- 
tained from recent calculations for several experimental 
high-speed airplanes have appeared much different from the 
conventional stability boundaries. Because of the com- 
plexity of the complete expression for 7?=0, it is difficult to 
determine the reasons for such unusual looking curves and 
the significance of the boundaries. From the simplified 
expressions, however, a complete analysis of the boundaries 
can be easily obtained. The 77= 0 boundaries of an experi- 
mental airplane are shown in figure 13 (a). In addition to 
the 77=0 boundaries, the 19=0 boundaries are also plotted 
in the figure. As mentioned at the outset of this report, 
77=0 is a neutral-oscillatory-stability boundary only if 
D is positive. The 77=0 boundaries on the hatched side of 
D — 0 are not therefore neutral-oscillatory-stability bound- 
aries. The boundaries for the same experimental airplane 
calculated from the simplified expressions are plotted in 
figure 13 (b). The 77, = 0 and D— 0 boundaries which are 
not neutral-oscillatory-stabilitv boundaries, as determined 
by the analysis presented in the section entitled “Validity of 
79= 0 and 7? i = 0 as Neutral-Oscillatory-Stability Bound- 
aries,” are shown as dash-line curves in the figure. In 
D== 0, the coefficient of Ci 0 becomes zero at 0^ = 0.056 
and, therefore, the D — 0 boundary approaches positive 
infinity in the second quadrant at 0,,^ = 0.056. As C„„ 
increases above 0.056, 79= 0 returns from negative infinity 
and appeal’s in the first quadrant. Similarly, 7?i = 0 ap- 
proaches negative infinity when C nfj is approximately equal 
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(a) R = = 

(h)/?i=0and 71=0. 

Figure IS.— Oscillatory stability boundaries for an experimental hifrh -speed airplane. 


to 0.25 since the coefficient of (\ g in 72, = 0, A,K X z— A^l 3 , 
is zero at this value of C„ g . Above C„ g of 0.25, /7, = 0 
returns from positive infinity and appears in the second 
quadrant. It is necessary to note that in figure 13 (a) the 
neutral-oscillatorv-stability boundary is one continuous 
curve; whereas in figure 13 (b) this boundary is composed of 
two sections, one section of 7A = 0 and the other section of 
71=0. The latter fact provides the important information 
that the period of the. oscillation which becomes unstable 
upon passing through the 71=0 boundary is comparatively 
longer than the period of the oscillation which becomes 
unstable upon passing through the 77, = 0 boundary. 


CONCLUSIONS 

The following conclusions were reached from a theoretical 
investigation of a simplified method for obtaining and 
analyzing the neutral-lateral-oscillatory-stability boundary: 

1 . A necessaiy condition for the lateral-neutral-oscillatory- 
stability boundary is that R=D{BC — AD) — B 2 E= 0, where 
A, B, C, D, and E are the coefficients of the lateral-stability 
equation. The expression for 72=0 is approximated by the 
expressions Z1=0 and R 1 =BC—AD=0. Criterions arc de- 
rived which, if satisfied, indicate that the approximate 
expressions satisfy the necessary and sufficient conditions 
for a neutral-oscillatory-stability boundary. 

2. If 71=0 and 7?i=0 approximate 72=0, the curve D—0 
represents the neutral-oscillatory-stability boundary for the 
oscillation which has a period comparatively longer than the 
period of the oscillation for which 72, = 0 is the boundary. 

3. In general, the results of the computations obtained 
from Ri = 0 and 71=0 show very good agreement with the 
results calculated by the exact expression for 72=0. Specifi- 
cally, the results of the investigation indicated: 

(a) An increase in the absolute value of the parameter A 2 , 
which is equal to (\ — 2C L K X 2 (where, ('„ is the yawing- 
moment coefficient due to rolling-angular-velocity factor, 
(\ is the trim lift coefficient, and K z is the radius of gyration 
in yaw), causes a destabilizing shift in the branches of 
72 = 0 approximated by 71=0 and R 1 = 0. 

(b) The branch of 72=0 approximated by D—0 mainly 
depends upon the parameter A., and the damping-in-roll 
derivative (\ . The product-of-inertia term K xz also appears 
in 71, but it has a negligible effect on the branch of 72=0 
approximated by 71=0. 

(c) An increase in the relative-density factor causes a 
destabilizing shift on the branch of 72=0 approximated by 
7?i = 0 but does not affect the branch of 72=0 approximated 
by 71=0. 

4. The neutral-oscillatory -stability boundaries computed 
from the simplified expressions show excellent agreement 
with the corresponding boundaries presented in NACA 
TN 1282. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 4, 194-8. 



APPENDIX A 

DERIVATION OF TEST FUNCTIONS A r AND Ad 


For a given value of C n ., selected from the probable range 
of C, lfi of the airplane for which the lateral-stability analysis 
is to be made, let 


Ei=riC^+r s 

d#, 

( C, ' / s)fi,=o = 

—r 2 

t>c l0 11 

— V 
T\ 

D=d 1 Ci 0 A-d 2 

dD j 

AC l0 

{^' l a)n = o 

AT 2=d 

di 

B 2 E=e,C l0 +e 2 

bB 2 E 

e ' 

{p l fd b'-e=(i 

<M I 

l b 


As shown in figure 14 the exact roots of R=R x D—B-E=0 
occur at the intersection of the straight line B 2 E with the 
parabola R { D. In the vicinity of the point Ci„—r, at which 
_K 1 = 0, the curve RiD is approximated well by a straight 
line tangent to the curve at Ci B —r, that is, 

RiD ~(^!) r 0 {Cl f>~ r) = (-'V/i+rA) (C l0 -r) 

If there is a root of R=RiD—B 2 E=Q near R x — 0 (that is, 
if B 2 E intersects R^D near the point r in fig. 14), then 

R= (.—r^di+ndi) (C tp —r) — eiC !f> — e 2 =0 



Figure 14.— Graphical representation of the roots of the i?=0 boundary. 


or 


Ad 


go.oi 


whichever is the larger. 

The expressions for r\, (l u and o for use in equations (Al) 
and (A2) are 


r i — 8m&(-Ti-^-xz A 2 A 3 ) 


Thus, the approximate deviation of a root of R = 0 from I?, =0 
is given by 


Ar = (C^—r) =~r 


r,e 2 


T\ {}\d 2 


—r 2 e i 

-r 2 cl i — e,) 


ei(e — r) 

r l d l (d—r)+e l 


(Al) 


where 


di — M5-A2 


e — 2ix 6 A l 2 C L C rtr 


A, = K A - 2 C„ r + 2 K x ' 2 K z 1 C Yfi + C lp KA 


If this deviation, Ar, is small, then R y = 0 is a close approxi- 
mation to one branch of 7? = 0. A suitable criterion for this 
approximation is 


A 2 =C„ p -2C L K/ 
ANKARA ~K X A 
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or 


Ar, 


^ 0.01 


whichever is the larger. 

In the case of D— 0, a similar analysis results in the test 
function 
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6i (e—d) 
f'ldi (r — d) +&i 


(A 2) 


If Ad is small, D= 0 may then be considered a close approxi- 
mation to the other branch of R— 0. A suitable criterion for 
this approximation is 
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